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Deriving exact results for Ising-like models from the cluster variation 
method 
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The cluster variation method (CVM) is an approximation technique which generalizes the mean field approxi- 
mation and has been widely applied in the last decades, mainly for finding accurate phase diagrams of Ising-like 
lattice models. Here we discuss in which cases the CVM can yield exact results, considering: (i) one-dimensional 
systems and strips (in which case the method reduces to the transfer matrix method), (ii) tree-like lattices and (iii) 
the so-called disorder points of euclidean lattice models with competitive interactions in more than one dimension. 



The cluster variation method (CVM) is a hi- 
erarchy of approximation techniques for discrete 
(Ising-like) classical lattice models, which has 
been invented by Kikuchi In its modern 

formulation |Q the CVM is based on the varia- 
tional principle of equilibrium statistical mechan- 
ics, which says that the free energy F of a model 
defined on the lattice A is given by 

F = min F[pa] = min Tr{pAH + pA In pa), (1) 

where H is the hamiltonian of the model, = 1 
for simplicity, and the density matrix pA must be 
properly normalized: Tt{pa) = 1. 

As a first step one usually introduces the cluster 
density matrices and the cluster entropies 



Pa 



Tr 



A\aPA 



-Tr(/9Q In Pa), 



(2) 



where a is a cluster of Ua sites and TrA\Q denotes 
a summation over all degrees of freedom except 
those belonging to the cluster a. One then in- 
troduces the cumulant expansion of the cluster 
entropies 
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(-ir°-"^5a, (3) 



in terms of which the variational free energy can 
be rewritten as 



F[pa] = Tr{pAH) - 
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(4) 



The above steps are all exact and the approx- 
imation defining the CVM comes in when one 
truncates the cumulant expansion of the entropy. 
The sum of the cumulants of the cluster entropies 
is restricted to a given set M of clusters, which in 
most cases can be thought of as a set of maximal 
clusters and all their subclusters. If the model 
under consideration has only short range inter- 
actions and the maximal clusters are sufficiently 
large the hamiltonian can be decomposed into a 
sum of cluster contributions and the approximate 
variational free energy takes the form 



F[{pa,aeM}] 
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[TiipaHa) - S 
[Tr{paHa) - aaS, 



(5) 



where the coefficients Cq. can be easily obtained 
from the set of linear equations 
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and the cluster density matrices must satisfy the 
following conditions which express normalization 



Trpa = 1, Vq G Af 
and compatibility 
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Having introduced an approximation it is 
worth asking whether there are special cases in 



2 



which it turns out to be exact. The simplest ex- 
ample is that of a system defined on a lattice A 
which can be regarded as the union of two clusters 
A = AUB, such that, denoting hy K = ACiB their 
intersection, there is no interaction between A' = 
A\K and B' = B\K. In this case the hamilto- 
nian has the general form H = HA{Q_A) + HB{g_B) 
and it is easy to check that the density matrix can 



be written as pA — 



PaPb 

PK 



, which in turn implies 



the decomposition — Sa + Sb ~ Sk for the 
entropy. The CVM approximation which one ob- 
tains with the set of clusters M = {A, B, K} leads 
to the same decomposition of the entropy (eq. ^ 
yields = as = 1, = —1) and is therefore 
exact. It can be verified that this argument can 
be easily generalized (to several clusters sharing a 
common intersection) and/or iterated (to a chain 
of clusters A, B,C, . . .). 

The above argument could be used to explain 
the well-known fact that CVM approximations 
are exact for Bethe and cactus lattices (that is, 
interior of Cayley and Husimi trees, respectively), 
which are made of links (respectively plaquettes), 
sharing common sites with no loops (respectively 
no loops larger than the elementary plaquette). 
However we shall leave apart tree-like lattices and 
turn our attention to euclidean ones, considering 
first one-dimensional systems (strips) and then 
the disorder points of higher dimensional models. 

Since it is known that the Bethe-Peierls ap- 
proximation (which is the lowest level CVM ap- 
proximation for a model with nearest-neighbour 
interactions only, obtained by taking M = 
{links, sites}) is exact for a one-dimensional 
chain, one might wonder whether there is a CVM 
approximation which is exact for a strip of finite 
width. The answer is affirmative and it is interest- 
ing to note that one recovers the transfer matrix 
formalism. Consider a strip of width N and (fi- 
nite, for the moment) length L and let the hamil- 
tonian contain only translation-invariant NN in- 
teractions, with open boundary conditions. In 
the above scheme, this is an example of a chain 
of intersecting clusters and we can guess that a 
CVM approximation with the N x 2 ladders as 
maximal clusters should be exact. Denoting by 
II such clusters and by I their x 1 intersections 



we set AI — {II, 1} (no other subclusters enter the 
cumulant expansion) and in the thermodynamic 
limit L — > oo, assuming that translational invari- 
ance is recovered, we get the variational principle 



lim — 

L — >oo Ij 



f = 

= min Tr {puHu 



piilnpii - pilnpi) . (9) 



Denoting by a and a' the two sets of degrees of 
freedom of the two I subclusters of a II cluster we 
can solve for pu and recover the transfer matrix 
formalism in the form 
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Inmax V pj/'(a)e-^"(^^^'Vi'^'(^') (10) 
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with the normalization constraint pi (a) 



1. 



It is interesting to note that the CVM comes 
with a natural fixed point algorithm |^ for find- 
ing the local miima of the free energy, which in 
this case reduces to the power method for finding 
the largest eigenvalue of the transfer matrix. 

The last (and perhaps the most interesting) 
case we want to consider is that of disorder points. 
As an example we consider the square lattice Ising 
model with competitive interactions, with hamil- 
tonian 

(i-j) 

-Ki ^ criCTjO-feCr/, (11) 
[ijkl] 

where if i > is the NN coupling, K2 < the 
next nearest neighbour coupling and the pla- 
quette coupling. It is known jj] that in the disor- 
dered phase of this model there is an integrable 
subspace given by 

w<..^^ e^^-i cosh(4X2) + e-2^2 

where the free energy density is given by 

/ = - In [exp(-X4) + exp(X4 - 2K2)] . (13) 

In this subspace the R matrix has an eigenvector 
which is a pure tensorial product and the eigen- 
vector of the transfer matrix corresponding to the 
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largest eigenvalue is also a pure tensorial prod- 
uct iQ], and hence the density matrix and the 
two-site correlations are factorized. Because of 
this factorization (and the corresponding decom- 
position of the entropy) one can expect that the 
model can be solved exactly by the CVM and in- 
deed this is the case. It is enough to choose M 
= {plaquettes and their subclusters} (of course 
larger maximal clusters work as well) to verify eq. 

and also to calculate the two-site correlation 
functions r(a;, y) = {(j{xo,yQ)a{xo +x,yQ + y)): 

T{x,y) = g\-\+\y\ 

cxp(-4ir2) - cosh(2Xi) 
^ sinh(2i^i) ^ ' 

as well as many-site correlation functions like the 
plaquette correlation function q ~ {(Tiajak'Ji)n 



e4if4 (1 _ e8^f2) + 4e2if2 [e^K^ + ^2X2) 



(15) 



Details, and generalizations to other models, will 
be reported elsewhere [|] . 
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